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Fluctuations in the electron system of a 
superconductor exposed to a photon flux 
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In a superconductor, in which electrons are paired, the density of unpaired electrons should 
become zero when approaching zero temperature. Therefore, radiation detectors based on 
breaking of pairs promise supreme sensitivity, which we demonstrate using an aluminium 
superconducting microwave resonator. Here we show that the resonator also enables the 
study of the response of the electron system of the superconductor to pair-breaking photons, 
microwave photons and varying temperatures. A large range in radiation power (at 1.54 THz) 
can be chosen by carefully filtering the radiation from a blackbody source. We identify two 
regimes. At high radiation power, fluctuations in the electron system caused by the random 
arrival rate of the photons are resolved, giving a straightforward measure of the optical 
efficiency (48 ±8%) and showing an unprecedented detector sensitivity. At low radiation 
power, fluctuations are dominated by excess quasiparticles, the number of which is measured 
through their recombination lifetime. 
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I n a superconductor well below its critical temperature, the 
I majority of the electrons is bound in a condensate of Cooper 
pairs. The further the superconductor is cooled down, the 
closer it gets to its ground state, where all the quasiparticles are 
condensed to pairs. Due to the low gap energy, the super- 
conductor is sensitive to disturbances from the environment to 
which it couples. In most experiments, this sensitivity is 
undesirable, but it is particularly suited for detection of radiation. 
The superconductor can interact with its environment due to 
either photons or phonons. Photons with an energy higher than 
the energy gap break up Cooper pairs into quasiparticles. The 
change in the number of quasiparticles and Cooper pairs changes 
the electrodynamic response of the superconductor, which can be 
measured using a microwave resonator 1 . Quasiparticles give 
rise to microwave losses and the Cooper pairs to a kinetic 
inductance 2 . In steady state, the number of quasiparticles 
fluctuates in time around a constant average value. A 
measurement of the spectrum of these fluctuations allows for a 
characterization of the quasiparticle system when exposed to pair- 
breaking photons, microwave photons or to changes in the bath 
temperature. The characteristic timescale of the fluctuations, the 
quasiparticle recombination time, is inversely proportional to 
the number of quasiparticles 3 , and is therefore a measure of this 
number. These fluctuation phenomena are a monitor of the 
superconducting state and reveal the physical mechanisms that 
are at the heart of pair breaking in a superconductor. 

We study these processes in a superconducting pair-breaking 
detector formed by a 50 nm thick Al film. The ideal pair-breaking 
detector can either count single photons while its sensitivity is 
limited by Fano noise, or is photon integrating and limited by 
photon noise, the noise from the photon source itself 4 . In both 
cases, a high optical efficiency is required. The principle of 
radiation detection due to pair breaking with superconducting 
microwave resonators was proposed in 2003 (refs 1,5). Since then, 



the most important source of excess noise , due to dielectric two 
level systems, has been thoroughly studied 6-9 . Here we report on 
an all- aluminium antenna- coupled microwave resonator detector 
(Fig. la) that is limited only by fluctuations in the electron system 
that are fundamentally connected to the physical process of pair 
breaking. 

Results 

Detector sensitivity. We use a blackbody with a variable tem- 
perature (3-25 K) and eight optical filters, which define an optical 
band around 1.54 THz, as shown schematically in Fig. lb. The 
radiation power can be changed from 3 x 10 _21 Wto7x 10 _ 13 
W. At powers ranging from 0.1-700 fW, the sensitivity is only 
limited by the random arrival rate of the photons, which is evi- 
dent through the measured power dependence of the noise 
equivalent power (NEP) as shown in Fig. lc. At lower radiation 
powers, we observe a power-independent NEP. This is consistent 
with generation -recombination noise due to the presence of 
excess quasiparticles 10 ' 11 . Excess quasiparticles are a general 
concern for superconducting devices 1 8 . In this case, they are 
generated by the microwave readout power 19,20 . As shown in 
Fig. lc, they impose a lower limit to the NEP of this detector of 
3.8 ± 0.6 x 10 ~ 19 WHz~ 1/2 , which is the lowest reported so far 
for this type of detectors. 

Design of the experiment. The detector is based on a lens 
antenna- coupled superconducting microwave resonator. The 
resonator is an open ended, half wave, coplanar waveguide 
resonator, where the central strip (with a width of 3 jim) is iso- 
lated from the ground plane. The resonators all have different 
lengths and therefore different resonant frequencies, enabling the 
readout of all resonators using a single coaxial line. Radiation is 
focused by a silicon lens onto an X-slot antenna 21 , optimized for 




Figure 1 | Overview of the experiment and main results, (a) A picture of one antenna-coupled microwave resonator. It consists of a half wavelength 
coplanar waveguide microwave resonator with two open ends, capacitively coupled to a microwave readout line. In the middle, it has an X-slot 
antenna to receive terahertz radiation. The ground plane and central strip of the resonator are respectively 100 and 50 nm thick layers of aluminium, 
(b) Schematic of the setup for measurements at various radiation powers. A blackbody with a variable temperature illuminates the lens 
antenna-coupled resonators through three stacks of filters, which define a passband around 1.54 THz. Because of the box-in-box configuration at 
100 mK and the coax cable filters, the device is well shielded from stray light, (c) The optical sensitivity is expressed in NEP as a function of radiation 
power at a frequency of 20 Hz. We observe two regimes: above 0.1 fW the NEP increases with y/P^, indicative of photon noise, whereas below 
0.1 fW the NEP saturates. The blue line is the photon noise limit as a function of power, with the optical efficiency (48%) taken into account. 
For the red dashed line, the generation-recombination noise limit due to excess quasiparticles is taken into account, based on a quasiparticle 
recombination time of 3 ms. The error bars (1 s.d.) are combined statistical uncertainties from the noise level and responsivity. 
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broad band detection from 1.4-2.8 THz. Radiation coupled to the 
antenna is launched as a travelling wave into the waveguide 21 , 
where it is absorbed by breaking Cooper pairs (the gap energy 
A= 188 ueV). The created quasiparticles, which can diffuse over 
several millimetres before they recombine, are confined to the 
central strip. The central strip layer is 50 nm thick, and the 
ground plane layer 100 nm. The thin central strip layer gives 
higher response and ensures that most of the radiation is 
absorbed in that central strip, due to its higher resistance (see 
Methods). The thick ground plane reduces antenna losses. An 
advantage of the geometry, shown in Fig. la, is that it can also be 
used at other radiation frequencies by only changing the antenna. 

The sample is cooled in a pulse tube precooled adiabatic 
demagnetization refrigerator. The sample stage is carefully 
shielded from stray light from the 3 K stage of the cooler, using 
a box-in-a-box concept with optical filters at each stage, as well as 
coax cable filters in the outer box 22 . The photon source is a 
blackbody with a variable temperature between 3 and 25 K 
(Methods). The system is schematically depicted in Fig. lb. Eight 
optical filters in series define an optical band-pass of 0. 1 THz 
centred around 1.54 THz. Three filter stacks are essential to 
eliminate filter heating. The filter transmission of the three filter 
stages is shown in Fig. 2a. The curves of spectral radiance for high 
and low blackbody temperature indicate a large tuning range in 
radiation power (P ra d)- In fact, P ra d can be varied between 3 zW 
(1 zW=10~ 21 W) and 1 pW (Supplementary Fig. SI). 
Practically, this experiment allows us to switch from a regime 
where the number of quasiparticles is fully determined by the 
radiation to a regime with a negligible number of optically created 
quasiparticles. We put a polarizing wire grid just before the 
detector to make sure the detector only receives radiation in the 
polarization direction of the antenna. 

Operation principle. The number of quasiparticles is measured 
through a measurement of the complex conductivity of the 
superconductor. The real part of the conductivity, o l9 is due to the 
quasiparticles and resistive. The imaginary part, cr 2 , is due to 
the kinetic inductance of the Cooper-pair condensate 2 . When the 
radiation power or the bath temperature is increased, more 



quasiparticles are generated, which changes both a 1 and cr 2 . The 
kinetic inductance i ncreases, wh ich leads to a lower resonant 
frequency f 0 = 1/2/^/ (L g + Lk)C, where / is the length of the 
resonator, L g the geometrical inductance, L k the kinetic 
inductance and C the capacitance of the line, all per unit 
length. The losses at microwave frequencies also increase, leading 
to a shallower resonance curve. Measurements of the resonance 
curves for various radiation powers are shown in Fig. 2b. In a 
practical detection scheme, one typically uses an amplitude, A, 
and a phase, 9, referred to the resonance circle in the complex 
plane 6 , as shown in Fig. 2c. The amplitude response originates 
from a change in resistance, whereas the phase changes due to the 
kinetic inductance. We have only used the amplitude response in 
this experiment, because the phase noise is too high for this 
resonator (Supplementary Note 1). 

The NEP is a convenient quantity to compare the spectra of 
quasiparticle fluctuations in different regimes, as shown in Fig. lc. 
The NEP of the resonator amplitude is experimentally deter- 
mined from a measurement of the noise spectrum (S A ) and the 
responsivity to radiation (dA/dP rSL d) and given by 

- i 



NEP(f) = y/sl(f) 



dA 
d?7 



rad 



l + ^Tl/Tqp) 2 , (1) 



with P rad the radiation power and / the modulation frequency. 
dA/dP r2L d is obtained experimentally by a linear fit to a 
measurement of A where P ra d is slowly varied around the power 
of interest (see Methods). The measured microwave response 
upon a change in P rad is shown in Fig. 2c. The last factor in 
equation (1) arises because the quasiparticle system cannot 
respond to fluctuations that are faster than the quasiparticle 
recombination time, Tqp. 

Photon-induced quasiparticle fluctuations. If the average 
number of quasiparticles is dominated by the absorbed optical 
photons, the number of quasiparticles fluctuates in time due to 
two contributions. One is fundamental to every power- integrating 
detector and due to the random arrival rate of the photons, which 
induces a random generation of quasiparticles. The power spec- 
tral density of fluctuations in the resonator amplitude due to this 
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Figure 2 | Response to radiation, (a) Filter transmission characteristics of the three stacks of optical filters in the setup (Fig. lb). The first and third 
set of filters have a low-pass, a band-pass and a high-pass filter. The second set (at the 100 mK box) has only a high- and a low-pass filter. In the 
bottom panel the total transmission of these eight filters is shown. We also show the normalized spectral radiance (Planck's law) at two blackbody 
temperatures, which demonstrate the large tunability in radiation power in this spectral range. Note that especially for low blackbody temperatures 
only a fraction 10 ~ 6 of the total power is in the spectral range of interest. The rejection of the rest of the power requires the eight consecutive filters, 
(b) The magnitude of the microwave transmission |S 2 i| 2 , measured as a function of frequency for various radiation powers as shown in the legend. 
At higher power, more quasiparticles are created, which give a higher resistance and inductance and therefore lead to a lower resonant frequency 
and a shallower dip. The dots show the resonant frequency at each power, (c) The resonance circle for a selection of radiation powers (legend), 
measured as a function of frequency (lines). The squares show the response upon a small change in the radiation power measured at constant 
frequency, the resonant frequency of each circle. In the last circle, we show how that response is translated into an amplitude, A, and a phase, 9. 
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photon noise is given by 



S ? A (f) = 2hFP md (l+t lovt B) 



1 + (27tfTqp) 2 ' 



(2) 



where the first term is the spectrum of the photon (power) 
fluctuations and the second term describes the resonator response 
upon a change in the radiation power. F is the frequency of the 
optical photons and h Planck's constant. The factor (1 -\-f] opt B) is 
the correction to Poissonian statistics due to photon bunching for 
a single mode, with f] opt the optical efficiency and B the mode 
occupation 23 . The photon bunching contribution is negligible 
for the here measured power range (Supplementary Fig. Sic). 
Equation (2) is valid as long as T qp »T res , which holds in this 
experiment since the response time of the resonator, given by 
x res = Q/nfo, is 6 jis. Q is the quality factor of the resonator. 

Because of the pair-breaking nature of the radiation 
absorption, a second noise mechanism arises due to random 
recombination of the quasiparticles that are generated by the 
photons. This is half the generation-recombination noise that 
arises in thermal equilibrium 24 , because generation noise is 
already contained in equation (2). The spectrum is given by 

2 



s R A (f) 



2ATqpTqp 



1 + (2tt/i 



qp/ 



dA 
dK 



qp 



(3) 



with IV qp the number of quasiparticles and dAldN qp the 
responsivity of A to a change in IVqp. Quasiparticle number 
fluctuations can be converted to power fluctuations through 
^pb^opt^rad = NqpA/T qp . rj opt is the optical efficiency, the efficiency 
with which power in front of the lens is absorbed in the detector. 
77 p b~0.57 is the pair-breaking efficiency 25 , the efficiency with 
which absorbed radiation power is converted into quasiparticles. 
For small changes in the quasiparticle number, dP TSL JdN qp = 
A/(T qp ^ pb ^ opt ) and therefore dA/dP rad = T qp r] ph r] opt /A • (dA/dN^). 
From equations (2) and (3), the relative contribution of 
photon noise compared with recombination noise is given by 
hF(\ + mB)r] ph r] opt / A =10 at all P rad , for F=1.54THz and 
rj opt = 0.5. 

The NEP due to photon noise and recombination noise 
(equations (1-3)), for/< l/(27iT qp ), is given by 



NEP, 



photon 



2P rad ^F(l+> ?opt B)+2AP rad /> ?pb 

'/opt 



(4) 



which is shown as the blue dashed line in Fig. lc. The three 
contributions in equation (4) are compared in Supplementary 
Fig. Sic. 

In thermal equilibrium IV qp is related to T qp through 



1 T 0 N 0 (hT c yv 



qp 



L qp 



2A 2 



(5) 



where N 0 is the single spin density of states at the Fermi level, 
V the volume, k B Boltzmann's constant, and t 0 the charac- 
teristic electron-phonon interaction time 3 . We take 
N 0 = 1.72 x 10 10 um~ 3 eV _1 and V=0.6x 10 3 um 3 , half the 
volume of the central strip of the resonator (see Methods). In this 
experiment equation (5) can be applied under the assumption 
that the response of the resonator amplitude is dominated by 
quasiparticle generation and recombination. We have shown 
earlier in a similar aluminium device that IV qp and T qp are related 
by equation (5) for thermally generated quasiparticles. This 
equation is also expected to hold in non-equilibrium conditions 
due to optical excitations 26 or due to microwave readout power 
dissipation at low bath temperatures 20 . 

The fluctuations in the resonator amplitude were measured as a 
function of radiation power at a constant bath temperature of 



120 mK (Methods). The resulting power spectral densities are 
shown in Fig. 3a for a selection of radiation powers and a 
constant microwave readout power (the power on the readout 
line) of — 88 dBm. We observe that the spectra show a roll-off, 
the frequency of which increases as a function of radiation power, 
due to the decreasing quasiparticle recombination time. The other 
phenomena in the noise spectrum at higher frequencies (the 
bump at 20 kHz, a second higher frequency roll- off and amplifier 
noise, see Supplementary Fig. S2) are understood and can be 
accounted for (see Supplementary Note 2). The quasiparticle 
recombination times from the roll- off in the spectra are shown as 
a function of radiation power in Fig. 3b. iV qp , calculated using 
equation (5) is shown on the right axis. Since ^pb^opt^rad — N qp A/ 
T qp and Ar q pOcl/T q p (equation (5)), T qp is expected to scale as 
T qp ocP rad / . A fit to the measured recombination time as a 
function of P rad results in TqpOcP rad 0 - 50±0 03 , which agrees very 
well with the expected behaviour. Within the measurement 
accuracy, the same coefficient is measured for other microwave 
readout powers. The quasiparticle recombination time saturates 
below about 0.1 fW. This saturation in the recombination time is 
commonly observed 27 , but not yet understood for all materials. 
Here the saturation is consistent with the presence of excess 
quasiparticles, which we have demonstrated in measurements on 
a very similar Al resonator 11 . 

Phonon-induced quasiparticle fluctuations. Excess quasi- 
particles give rise to quasiparticle number fluctuations 11 ' 19 . 
To verify that the spectra in the saturation regime show these 
fluctuations, we change the number of quasiparticles by varying 
the number of phonons (the bath temperature) at the same 
microwave power. The amplitude spectrum is shown for bath 
temperatures ranging from 90-255 mK in Fig. 3c. The blackbody 
temperature is kept at 3.2 K, so there are less than 100 
quasiparticles due to the radiation power in the sensitive 
volume. The amplitude spectrum due to quasiparticle fluctu- 
ations can be described as 10,11 S^ R (f) — 2S^(f), because here 
both generation and recombination are considered. The noise 
level, which is proportional to iV qp Tqp, is expected to be constant 
as a function of temperature (see equation (5)), which we indeed 
observe in Fig. 3c. We assume here that dA/dN qp is constant for 
this temperature range 26 . The quasiparticle recombination time 
extracted from these spectra is plotted as a function of 
temperature in Fig. 3d. We observe the same saturation level in 
the recombination time as in Fig. 3b where the radiation power 
was changed. t 0 = 303 ± 14 ns is obtained from a fit to the 
measured i qp as a function of temperature 3 ' 11 . t 0 is slightly 
different from earlier results (458 ns 11 ), which could be due to the 
higher resistivity and T c of the Al 3 . 

We have now verified that the noise spectra in the regime 
of low radiation power (below 0.1 fW) are consistent with 
quasiparticle number fluctuations. The optical NEP due to 
quasiparticle number fluctuations is given by 



NEP GR 




(6) 



Using this equation, the saturation of the NEP in Fig. lc is 
consistent with the measured saturation in the quasiparticle 
recombination time (Fig. 3b) and the number of quasiparticles 
inferred from that measurement. 

If we return to the photon -induced fluctuations in Fig 3 a, we 
observe that the noise level becomes constant at the highest 
radiation powers as well. This constant level is expected when 
equations (2) and (3) are rewritten in terms of NqpT q p. The noise 
level is higher than in Fig. 3c, because there is both photon noise 
and recombination noise. 
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Figure 3 | Quasiparticle fluctuations, (a) Power spectral density of the resonator amplitude as a function of frequency for different radiation powers at a 
constant bath temperature of 120 mK and a constant microwave readout power of -88dBm. Lorentzian fits to the spectra at the lowest and highest 
temperatures (dashed lines) show how the quasiparticle recombination time can be extracted from the spectra. A noise floor due to amplifier noise is 
added to the fitted roll-off. (b) The quasiparticle recombination time as a function of radiation power obtained from the roll-off frequency in the 
measured spectra. The error bars denote statistical uncertainties from the fitting procedure (1 s.d.). The line is a power law fit to the last five points 
(where i qp does not saturate): Tq P ocP r ~ d 0 - 50±0 03 . The right axis shows the number of quasiparticles corresponding to the measured recombination time. 

(c) Power spectral density of the resonator amplitude as a function of frequency for different bath temperatures at a microwave readout power of 
-88dBm. As expected, the level of the spectrum stays constant and the roll-off frequency increases with increasing temperature, corresponding to a 
decreasing recombination time. At the highest two temperatures, the spectral level starts to rise, because the amplifier noise starts to dominate. 

(d) Quasiparticle recombination time as a function of temperature as extracted from the spectra. The error bars are 1 s.d. The solid line is the theoretical 
expectation for the recombination time from ref. 3. The right axis shows the number of quasiparticles corresponding to the measured recombination time. 



Noise equivalent power. The measured NEP, obtained by using 
equation (1) together with the measured S A , dA/dP rad and i qp , 
is shown for various radiation powers in Fig. 4. The NEP 
measurement was done at a range of microwave readout powers. 
The results shown in Figs lc and 4 are at the readout power with 
the minimum NEP for that radiation power. The measured dAI 
dP rad and S A that constitute the NEP are shown in the 
Supplementary Fig. S3. Note that the levels of these S A are dif- 
ferent from those in Fig. 3a because of the strong dependence of 
dA/dP r:i d on the microwave readout power. 

The measured optical NEP at 20 Hz is shown as a function of 
radiation power in Fig. lc as our main result. At radiation powers 
of 0.1 fW and higher, the NEP scales with \/P rSL d> as expected 
from the photon noise limit given by equation (4). In this regime, 
the optical efficiency is obtained by fitting equation (4) to the 
measured NEP. The result is shown as the blue line in Fig. lc, 
which gives rj opt = 0.48 ± 0.08 for a single polarization, consistent 
with electromagnetic simulations of the antenna (Methods). 

Below 0.1 fW, the NEP saturates at 3.8 ± 0.6 x 10" 19 WHz~ 1/2 . 
We have seen that generation-recombination noise due to excess 



quasiparticles dominates the noise spectra in this regime. From 
the measured recombination time (3 ms, see Fig. 5b), we calculate 
iVqp using equation (5). The sum of equations (4) and (6) is 
shown as the red dashed line in Fig. lc and gives a good account 
of the measured NEP. The limit of 3.8 x 10" 19 WHz~ 1/2 is in 
good agreement with predictions based on dark experiments 19 . 

Excess quasiparticles due to the readout power. From 
recent dark experiments 19 (without blackbody source) and 
simulations 20 , we expect that the excess quasiparticles at low 
radiation power are due to microwave readout power dissipation. 
In an optical experiment, excess quasiparticles should lead to a 
decrease in the optical response. The measured responsivity is 
shown in Fig. 5a as a function of P rea d- Indeed, the responsivity 
decreases with increasing P rea d> consistent with excess quasi- 
particles due to microwave absorption. A power law fit to the 
responsivity versus readout power results in dA/P rad ocP r ~ ad -4±01 . 
Figure 5b shows i qp , extracted from the noise spectra, as a 
function of P rea d- T qP increases when P read decreases, which is also 
consistent with quasiparticle generation by the micro - 
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Figure 4 | Noise equivalent power. Measured optical NEP in the resonator 
amplitude as a function of frequency for different radiation powers. The 
corner frequency of each spectrum corresponds to the quasiparticle 
recombination time as plotted in Fig. 3b. The stars indicate the NEP at the 
reference frequency of 20 Hz, which is chosen well within the quasiparticle 
roll-off. These are the NEP values shown in Fig. 1c. The measurements 
shown are taken at the readout power that gives the minimum NEP for that 
radiation power. 

wave readout signal 19 ' 20 . A fit to the measured data gives 

~ ^p- 0.2 ±0.05 
%>°^ re ad 

Microwave absorption leads to a strong non-equilibrium 
quasiparticle distribution in aluminium microwave resonators, 
which was recently shown in calculations 20 and experimentally 28 . 
The complex conductivity is shown to be more sensitive to the 
details of the non-equilibrium quasiparticle distribution than i qp . 
Therefore, the different exponents for dA/dP rad and Tqp in Fig. 5 
are not surprising. An extension of the work in ref. 20, including 
pair-breaking radiation, would be needed to understand the data 
in Fig. 5. If we make two assumptions, we can get a rough 
estimate of how the observables in Fig. 5 scale with P rea d- We 
assume that the absorbed microwave power is proportional to 
P rea( i and that ^ rea d> the efficiency with which the readout power is 
converted to quasiparticles, is constant. In that case, the optical 
responsivity is expected to scale with readout power as 
dA/P md ocP rea ^ -5 in the regime where the readout power dominates 
the number of quasiparticles, and ^qp^P^ 5 ( see Supplementary 
Note 3 for a derivation). Although the scaling of dA/dP rad is in 
reasonable agreement with the measurements, the approximation 
does not cover the different exponents for responsivity and t^. 
A more detailed theoretical treatment of this problem is required. 

Discussion 

At higher radiation powers, where the noise spectrum is domi- 
nated by photon noise, the optical responsivity also changes with 
readout power. In this regime (P ra d>l fW) however, the 
measured photon noise NEP stays the same, as expected from 
equation (4) (Supplementary Fig. S4). Therefore, when photon 
noise dominates the noise spectrum, one can safely use high 
readout powers to suppress amplifier noise (Supplementary Note 
4). 

At the lowest readout power where i qp was determined, 
^Vead = — 102 dBm, the quasiparticle recombination time is 3.5 ms, 




1.5 



-110 -105 -100 -95 -90 
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Figure 5 | Readout power dependence, (a) The optical responsivity, as 
obtained from a fit to the measured amplitude response, as a function of 
microwave readout power. The error bars denote statistical uncertainties 
from the fitting procedure (1 s.d.). The line is a power law fit to the data: 
oc ^read^ ^ ^ ■ Cb) The quasiparticle recombination time as obtained 
from the roll-off in the noise spectra as a function of microwave readout 
power. The error bars are 1 s.d. Below —102 dBm, the roll-off due to the 
recombination time is not visible anymore. The line is a fit to the data with a 
power law coefficient of Tqp°cP r ead' 2±0 ' 1 - AH data are measured at a 
radiation power of 3 zW and a bath temperature of 120 mK. 

which corresponds to a quasiparticle density n qp = 24um -3 . 
This density is still high in comparison with the lowest 
reported values for superconducting qubits and Cooper-pair 
transistors 17 ' 18 (less than 0.1 urn - 3 ), but inherent to the relatively 
high microwave powers we need in this type of experiments. The 
measured limit in optical NEP due to excess quasiparticles is 
comparable with the lowest observed optical NEP in other 
detectors for similar wavelengths 29-31 . 

A reduction in ATqp is possible by using a parametric amplifier 
with high bandwidth and dynamic range 32 . This allows a 
reduction of the readout power by about a factor 10. In the 
current design, however, the detector would become too slow for 
practical use at low readout power due to the long recombination 
time. The most feasible route towards lower NEP with 
aluminium, the most reliable material so far, is to choose 
geometries in which the active volume is dramatically reduced, 
which could also be the route towards single photon counting at 
terahertz frequencies. 

Methods 

Sample design. A layer of aluminium with a thickness of 100 nm is sputtered onto 
a sapphire substrate and serves as the ground plane for the microwave resonators. 
The microwave resonator is a coplanar waveguide resonator with a central strip 
width of 3 um and slit widths of 1.5 um. The central strip of the resonator is made 
of a second layer of 50 nm thick Al. The critical temperature of the 50 nm layer is 
measured to be T c = 1.24 K, from which the energy gap A = l.76k B T c =188 ueV. 
From the normal state resistivity, (p = 2.2 |if2cm for the central strip and 
0.28 |iQcm for the ground plane) the skin depth for radiation at 1.54 THz is 60 nm 
in the central strip and 21 nm in the ground plane. The X-slot antenna would be 
ineffective for a layer thinner than the skin depth, therefore the ground plane layer 
is 100 nm thick. The microwave sheet resistance of the central line is 0.46 Q and 
that of the ground plane 0.13 fl Taking this effect into account, together with the 
participation ratios of the central line and the ground plane 33 , we estimate that 90% 
of the radiation is absorbed in the central line. 
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The current distribution along the length of the resonator peaks at the antenna 
and decreases as sin(x) to zero at the open ends. Therefore, the responsivity 
changes with sin 2 (x). Since the diffusion length within a typical quasiparticle 
recombination lifetime of 2 ms is more than half the resonator length, optically 
created quasiparticles can move into the non-responsive regime. Therefore, for 
calculating the number of quasiparticles in the sensitive volume, we take half the 
central strip volume, V=0.6 x 10 3 um 3 . 

Noise measurement. The signal from the microwave generator is first attenuated, 
sent through the sample and amplified with a HEMT amplifier at 4 K and with a 
room temperature amplifier. The output is mixed with the original signal using an 
IQ mixer, the output of which can be sampled at a maximum frequency of 2 MHz. 
The spectrum of fluctuations in the resonator amplitude is measured by recording 
the resonator amplitude as a function of time and computing the power spectral 
density. Peaks in the time domain stream that occur due to high energy impacts are 
filtered out before the spectrum is computed, as described in ref. 11. We use the 
amplitude direction because fluctuations in the phase direction are dominated 
by two level system noise in the dielectrics surrounding the resonator 6 
(Supplementary Fig. S5 and Supplementary Note 1). 

Optical system and radiation power. The photon source used in the experiment 
is a blackbody that is formed by a 40 mm diameter copper cone, coated with 
carbon-loaded epoxy (EPOTEK 920 1LB part A, with 3% by weight carbon black 
and 3% by weight EPOTEK 920 1LB part B) that is covered with 1 mm SiC grains. 
The temperature of the blackbody is varied in this experiment from 3-25 K. 
There are three metal-mesh filter stacks (QMC Instruments, Cardiff), the 
characteristics of which are given in Fig. 2a. The measured transmission of the 
whole filterstack as a function of optical frequency, Tr{F), is also shown in Fig. 2a. 

Since there is no aperture limitation in between the blackbody and the detector, 
the optical throughput is assumed to be (c/F) 2 , with c the speed of light. The total 
radiation power that arrives in front of the lens of the detector can now be 
calculated by numerically integrating Planck's law over the throughput and the 
measured filter characteristic at each blackbody temperature T BB . The radiation 
power is here given for one polarization. 



Prad(T BB ) : 



JO 



Tr(F)hF dF 
exp(hF/k B T m ) - I 



(7) 



The optical window around 1.54 THz, together with the blackbody temperature 
range of 3-25 K gives a large tuning range in radiation power, as shown in 
Supplementary Fig. Sib. With the present device, we can verify the radiation power 
down to 0.1 fW using the measured quasiparticle recombination time (Fig. 5b). 
The trend of increasing optical responsivity and recombination time with decreasing 
microwave power, without any sign of saturation, suggests that the optical system is 
well characterized down to even lower radiation powers. The excess quasiparticles in 
the present device limit us to verify that. 

Right in front of the detector, after the last optical filter, we place a polarizer to 
select the polarization for which the antenna is designed. The polarizer consists of a 
copper wire grid on top of a 1.5-um thick Mylar film. The grid lines are 10 um wide 
and the spacing between the lines is 20 um. 

Responsivity fit. The responsivity dAldP r2L & is obtained from a linear fit to the 
response A measured as a function of P rad . The blackbody temperature is swept 
slowly around the P rad of interest. Each measured point is integrated for 500 ms. 
At very low radiation power (P rad <0.3 aW) there is no response (due to the excess 
quasiparticles). Therefore, the responsivity at 3 zW is estimated by a fit of the 
response in the range 0.3-9 aW, for example. Note that this reflects the saturation 
in the sensitivity of the detector. At the highest power P rad = 724 fW, the fit range is 
705-743 fW. The fit range is chosen by hand, because it differs widely from the 
saturation regime to the regime in which the number of quasiparticles is dominated 
by the absorbed P rad . The fit range is chosen small enough to stay in the linear 
response regime and large enough to limit the standard deviation of the fit. 

Lens antenna efficiency. The radiation power is focused by an elliptical silicon 
lens of 2 mm in diameter onto the antenna that is in the second focus of the lens 34 . 
The major and minor axes of the ellipse are 1.037 and 0.992 mm, respectively. The 
lens has an antireflection coating of 130 um of Parylene C, which is not optimized 
for 1.54 THz. The antenna is an in-line X-slot antenna, designed to receive 
radiation in a broad band around 1.54 THz as described in ref. 21. To obtain the 
optical efficiency, a simulation in CST Microwave Studio is performed of the whole 
structure: the lens with the coating, the antenna and a piece of coplanar waveguide 
transmission line. The aperture is chosen to be 30 degrees, the angle from which 
the detector can see the blackbody (single- side angle). The optical efficiency is 
shown as a function of frequency in Supplementary Fig. S6. The total efficiency is 
the multiplication of the front-to-back ratio, the spill over losses, the efficiency of 
an impedance mismatch between the antenna and the CPW line and the reflection 
losses at the anti-reflection-coated lens surface. All together, we expect an efficiency 
of 48% for one polarization in the filter transmission band, which is in good 
agreement with the measured optical efficiency of 48 ± 8%. The antenna was not 
designed to have a perfect optical efficiency, but to have a large bandwidth. The 



agreement of the measured optical efficiency with the CST- simulation shows that 
the optical system is understood. 
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